Mixtures of Bose gases confined in concentrically coupled annular traps 
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A two-component Bose-Einstein condensate confined in an axially-symmetric potential with 
two local minima, resembling two concentric annular traps, is investigated. The system shows 
a number of phase transitions that result from the competition between phase coexistence, and 
radial/azimuthal phase separation. The ground-state phase diagram, as well as the rotational prop- 
erties, including the (meta)stability of currents in this system, are analysed. 
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I. INTRODUCTION 

The field of cold atoms has expanded dramatically over 
the last 15 years. It has now reached a stage where ex- 
perimentalists are capable of designing the form of the 
confining potential. Going to extreme aspect ratios, con- 
ditions of quasi-one- and quasi-two-dimensional behavior 
have been achieved. In other experiments, it has also 
become possible to design toroidal trapping potentials 
0, & IS EL H, in which persistent currents have been 
observed 

Recent theoretical studies have examined Bose- 
Einstein condensates in one-dimensional annular traps. 
For example, quantum-tunncling-related effects in verti- 
cally [(| and concentrically 0] coupled double-ring traps 
were investigated. Also, the rotational properties of a 
mixture of two distinguishable Bose gases that are con- 
fined in a single ring were addressed Q . One of the basic 
points of the above studies is the fact that the ability to 
design traps, control and manipulate the atoms with a 
very high accuracy, may allow the investigation of novel 
quantum phenomena, like quantum phase transitions, for 
example. 

In the present work we consider a mixture of two dis- 
tinguishable Bose gases d, [l(| which interact via an 
effectively-repulsive contact potential, and are confined 
in a two-dimensional concentric double-ring-like trap, as 
shown in Fig.[TJ Using the mean-field approximation, we 
investigate two main questions: First, we identify the 
various phases in the ground state of the system, varying 
the interaction strength between the atoms. In the trap- 
ping potential that we consider, we observe that the two 
gases separate radially via discontinuous transitions; in 
this case, each gas resides in one of the two minima of the 
trapping potential, preserving the circular symmetry of 
the trapping potential. We also observe the expected az- 
imuthal (and continuous) phase separation between the 
two gases in each potential minimum [TJ [l2|, EH • A sim- 
ilar effect has also been studied in the case of a single 
ring @; see also p,Eip. 

The second main question that we examine are the ro- 
tational properties of this system, including its response 
to some rotational frequency of the trap Q, as well as the 
stability of the persistent currents for variable couplings, 



and variable relative populations of the two components. 
The expectation value of the angular momentum of the 
system as a function of £1 shows an interesting structure, 
reflecting the various phase transitions that take place 
with increasing fl. 

Regarding the (meta)stability of the currents, it is re- 
markable that for equal populations between the two 
components the vast majority of the coupling strengths 
that we have examined yield metastable states, except 
for a very small range where all the coupling strengths 
are exactly or nearly equal. 

It is worth mentioning that analogous single and con- 
centric ring geometries have been addressed in semicon- 
ductor heterostructures, both theoretically and experi- 
mentally, see e.g. 14 1, and also [lj| for reviews on the 
subject. In these systems, the applied external mag- 
netic field plays the same role as the trap rotation in 
the present problem and allows the investigation of, e.g., 
electron localization effects and persistent electron cur- 
rents in field- free regions. 

In what follows we first describe in Sec. II our model. In 
Sec. Ill we present the results for the ground state of the 
system, identifying the states where the species coexist, 
or separate, either radially or azimuthally. In Sec. IV we 
examine the rotational properties for a fixed rotational 
frequency of the trap, and the (meta)stability of the per- 
sistent currents. We study the stability as a function of 
the coupling between the atoms, as well as of the ratio 
of the populations of the two components. Finally, in 
Sec. V we present a summary and our conclusions. 



II. MODEL AND METHOD 

We consider two distinguishable kinds of bosonic 
atoms, labelled as A and B, which are trapped in a two- 
dimensional potential of the form 



V( P ) 



1 



min <^ - Mut (p - Rif , - Mu? Q (p - R ) 



(1) 

where p is the usual radial coordinate in cylindrical co- 
ordinates and M is the atom mass, assumed to be equal 
for the two components. The two (overlapping) parabo- 
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FIG. 1: (Color online) The confining potential V of Eq. JTJ, 
where Ri = 2oo, Ro = 4ao, and ui /cui = 5/4. 

lae in V(p) with frequencies and w are centered at 
the positions with p = Ri and p = R a , giving rise to the 
potential plotted in Fig.Q] . 

In our calculations we consider Ri = 2ao and R a = 4ao, 
where a = \h/ (Muj)] 1 / 2 is the oscillator length corre- 
sponding to uj = Wj/4, and finally L0 o /uJi = 5/4. In the 
outer ring the potential is more tight, oj > uii, in order 
to compensate for the fact that R Q > Ri, i.e., to make the 
product of the "width" of each annulus times the radius 
of each annulus to be comparable to each other. (This 
is typically also the case in the studies on electrons in 
quantum rings, in semiconductor heterostructures that 
we mentioned above.) 

To simplify the discussion, we also assume that there 
is a very tight trapping potential along the z axis (omit- 
ted in the potential above), which completely freezes out 
the degrees of freedom of the gases along this direction. 
With this assumption, our problem becomes effectively 
two-dimensional, with the tight dimension entering only 
implicitly through the parameters «,j in the Hamiltonian 
of Eq. ^j. With the usual assumption of a contact inter- 
atomic potential, the Hamiltonian becomes 

N A n 2 N B h 2 

t=l j=l 

N A N b 

N A ,N B 

+uab S ( T i- T i)-Q) 

i=l,j=l 

Here uy = (AttTi 2 aij / M) J |i/>o(.z)| 4 dz, with ipo(z) being 
the state of lowest energy of the potential along the z 
axis and ay being the s— wave scattering lengths for zero- 
energy elastic atom-atom collisions. The coupled Gross- 
Pitaevskii-like equations for the order parameters of the 
two components 0A and 4>b, resulting from the above 
Hamiltonian, are 

fi 2 V 2 

2M~^ A + V<kP ^ a + 9aa\4>a\ 2 4>a + ffAs|0s| 2 0A 

= (J-a4>a, 

h 2 v 2 

^TT^ B + V (P) ( I ) B + gBB^B^^B + gAB\4>A\ 2 4>B 



= p B <t>B- (3) 

In the above equations f \4>A\ 2 dr = Na/Nb, and 
/ \4>B\ 2 dr = 1. Also, gAA = N b uaa, 9bb = N b u B b, 
9ab = NgUABi an d Mi are the chemical potentials. 
In what follows, we consider repulsive interactions only, 
gu > 0, gij > 0, and also assume that gAA — gBB = .9- 

The method that we adopt to solve Eqs. ([3]) is a fourth- 
order split-step Fourier method within an imaginary-time 
propagation approach [l7|. We start with a reasonable 
initial state for the two components and propagate it 
in imaginary time, making sure that we proceed a suffi- 
ciently large number of time steps, which guarantee that 
we have reached a steady state. 



III. PHASE DIAGRAM OF THE GROUND 
STATE 

We start with the ground state of the system. There 
are three energy scales in the problem, namely the single- 
particle energy that is set by the trap, the intra-atomic 
interaction energy, and the inter-atomic interaction en- 
ergy. For weak interactions, the energy is dominated by 
the single-particle term and the ground state is deter- 
mined by the minimization of this term. On the other 
hand, for strong interactions, for a given ratio of the two 
populations Na/Nb, the actual symmetry of the ground 
state results from the competition between the intra-and 
inter-species coupling strengths; the former favors the 
maximum possible spread of the two gases within the sys- 
tem, whereas the latter favors the minimization of their 
spatial overlap. 

The most pronounced difference of this problem as 
compared to the case where there is only one potential 
minimum - i.e., a single annulus - is the existence of a 
phase where each component resides in only one of the 
two potential minima, thus separating radially. In addi- 
tion, we have also observed the expected azimuthal phase 
separation when both species occupy the same potential 
minimum @, M, E3, Ell- 
in Fig. [5] we illustrate the phase diagram showing the 
symmetry of the ground-state density distribution of the 
two gases as gAB and g are varied, for Na = Nb- Since 
gAA = gBB and Na — Nb, the results are symmetric 
when the two components are interchanged. The axes 
in the phase diagram of Fig. [2] may also be considered to 
represent the values of the scattering lengths clab and 
claa — o-bb (scaled appropriately). 

We have found three different phases: (i) coexistence 
of the two species (squares, red color), (ii) azimuthal sep- 
aration (triangles, blue color), and (hi) radial separation 
(circles, green color), see Fig. [3] The difference between 
solid and empty symbols in Fig. prefers to the stability of 
the persistent currents and is explained in the following 
section. 

As seen in Fig. [21 the phase boundary between phase 
coexistence and separation of the two components is, to a 



■ ■■■■■■■■ EnA i 

■ ■■■■■■■ eoA A 

■ ■■■■■■ EnAA ▲ 

■ ■■■■■ ehA A A A 

■ ■ ■ ■ ■ tttA a AAA 

■ ■ ■ ■ EnAA ▲ ▲ ▲ ▲ 

■ ■ ■ EnAA ▲ ▲ ▲ ▲ ▲ 

■ ■ EnAA ▲▲▲▲▲▲ 

■ EnAA A A A A A # 



■ ehA ▲ 
■ ehA A ▲ 

EnAA A A 

4 fb A ft 



A 



g AB (^M) 



FIG. 2: (Color online) Phase diagram in the (jAB-versus-g 
plane, which shows the symmetry of the ground-state den- 
sity distribution of the two components. Squares, triangles 
and circles correspond to phase coexistence, azimuthal phase 
separation, and radial phase separation, respectively. The 
solid (empty) symbols denote points where the currents are 
metastable (unstable). 

very good approximation, a straight line given by (setting 
h = M = uj = 1) 



9AB 



(4) 



as we have found by fitting numerically our data. 

Certain limiting cases in the phase diagram of Fig. 2 
may be analyzed and understood easily. In the case 
(jab = and fjAA = 9bb, the two gases do not interact 
with each other. In order to minimize their energy, they 
distribute homogeneously along the rings and coexist. In 
the other limiting case qaa = 9bb = 0, for gAB smaller 
than w 2, the two components also coexist; however, 
for larger values than 2, they separate azimuthally. Al- 
though in this phase the kinetic energy increases, the in- 
teraction energy is lowered due to the repulsion between 
the two species, and the azimuthal symmetry-breaking 
persists with increasing gAB- 

It is also instructive to understand the internal struc- 
ture of the phase diagram. As one moves vertically, i.e., 
for a fixed value of gAB (being sufficiently large, such that 
the components separate azimuthally), for small enough 
g the two components occupy mainly the inner ring since 
the repulsion between the particles cannot compensate 
for the stronger confinement of the outer ring. As g in- 
creases, the outer ring becomes progressively more occu- 
pied. When the inter-component repulsion gAB becomes 
large enough, the gases minimize their energy by sepa- 
rating radially. However, if g increases further, the inner 
ring becomes too small to host one of the species entirely. 
Azimuthal phase separation takes place again, now with 
both gases being largely spread within the whole system. 
Eventually, for even larger values of g, the dominant term 
in the energy is the intra-atomic interaction, and thus the 
two components coexist. 




FIG. 3: (Color online) Ground-state density distributions 
n A (p,0) = \(j> A (p,0)\ 2 and n B {p,6) = \4> B (p,0)\ 2 of the two 
components, for gAB = 55 and (a) g — 5, (b) g = 15, (c) 
g — 35, and (d) g = 60, corresponding to azimuthal phase 
separation, radial phase separation, azimuthal phase separa- 
tion, and phase coexistence, respectively. The peaks in the 
densities correspond to the minima of the inner and the outer 
parabolae of the confining potential. 



As anticipated before, we illustrate this effect in Fig. [31 
where we plot the densities of the cases with gAB = 55, 
and (a) g = 5, (b) g = 15, (c) g = 35, and (d) g — 
60, corresponding to azimuthal phase separation, radial 
phase separation, azimuthal phase separation, and phase 
coexistence, respectively. 

It is also of interest to investigate the nature of the 
phase transitions occurring in the system. As one crosses 
the boundary from coexistence to azimuthal phase sep- 
aration, the two components decrease continuously their 
overlap, developing sharper profiles as the repulsion in- 
creases. This transition is thus continuous (second or- 
der) , as it is also the case in purely one-dimensional single 
rings @. In the corresponding energy surface, the min- 
imum (which determines the ground state) moves con- 
tinuously as one crosses the phase boundary. On the 
contrary, the transitions involving radial separation are 
discontinuous (first order), indicating that two local min- 
ima in the energy surface compete and that the system 
jumps abruptly from the one state to the other. 



IV. ROTATIONAL PROPERTIES 

A. Lowest state of the system for a fixed angular 
frequency of rotation 

The phase transitions that we described above also 
have a clear influence on the rotational properties of the 
system. We start by examining the response of the sys- 
tem to some finite rotational frequency Q of the trap. In 
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the following we determine the total angular momentum 
per particle I — (L A + L B )/(N A + N B ) as a function of 
fi, I = /(fi). Following the usual procedure, we minimize 
the energy of the system in the rotating frame, i.e., we 
minimize E/N — /fi, where E is the total energy. The re- 
sult of this calculation is shown in Fig. [4) where we have 
set N A /N B = 1.5, g = 10, and g AB = 20. The typical 
values of fi are rather small because we have scaled it 
with u>i. However, at least in the case of a purely one- 
dimensional ring potential, the scale for the typical fi is 
on the order of the frequency fi/ (2MB 2 ) of the corre- 
sponding kinetic energy Ti 2 /(2MR 2 ) where R is the 
radius of the ring. In atomic units, h = M = uj = 1, 

= 4, while h/(2MR 2 ) ~ 1/50 (for R, say, equal to 5). 
This introduces a factor 10~ 2 . 

It is also instructive to comment on the behavior of the 
gas as fi increases and see the connection of the function 
I = /(fi) with the density distribution of the two compo- 
nents. For < fl/u>i < 0.005 the two species are sepa- 
rated radially and the total angular momentum is zero. 
For fi/wj 1=3 0.005, the density of the two components 
breaks its azimuthal symmetry discontinuously, and the 
angular momentum jumps abruptly to a finite value. Be- 
yond this value of fi ju)^ the angular momentum increases 
linearly with fi and is carried by both components, which 
undergo solid-body rotation. This is an expected result 
due to the azimuthal symmetry-breaking of the density, 
as we have confirmed by studying the phases of the two 
order parameters. When fi/cj, w 0.035, a new plateau 
appears, which corresponds to an angular momentum per 
particle equal to two, and the two species separate radi- 
ally. As one can see from the plot, there is a sequence 
between plateaus at integer values of / and straight lines 
with a positive slope, which are separated by abrupt 
jumps. This sequence persists up to fi/w^ ~ 0.085. Be- 
yond this value of fi there is no longer radial separation 
of the two components due to the large centrifugal force, 
which forces both gases to occupy the two potential min- 
ima and therefore to separate azimuthally. 

Similar discontinuous transitions in the function I = 
/(fi) occur in a single-component weakly-interacting 
Bose-Einstein condensate that rotates in a harmonic 
trap [l9l |. and are associated with discontinuous tran- 
sitions between phases of different symmetries of the 
single-particle density distribution of the gas. In the 
present problem, the corresponding different symmetries 
are the ones where the components separate radially, or 
azimuthally. 

Furthermore, the jumps in the I = /(fi) plot are consis- 
tent with the fact that for the parameters considered, the 
system supports persistent currents. In other words, had 
this function been continuous, then metastability would 
have not been possible (for a discussion of this effect we 
refer to Ref. [18||). 

From the above observations we see that the general 
picture that emerges by considering a finite rotation of 
the trap resembles the one where the couplings are var- 
ied, with a series of phase transitions between radial and 




a/co. 

in 

FIG. 4: (Color online) The total angular momentum per par- 
ticle, (La + Lb)/(Na + Nb), as a function of the angular 
frequency of rotation of the trap fi, which results from the 
minimization of the energy in the rotating frame. 

azimuthal separation. 

B. Metastable currents 

Another interesting question is the possible existence 
of (meta)stable currents. In physical terms, we inves- 
tigate the energetic (meta) stability of current-carrying 
states. In the case that there is an energy barrier that 
separates a current-carrying state from the ground state, 
in the presence of some dissipative mechanism (such as a 
thermal cloud, for example) such a state does not decay, 
and therefore the system supports persistent currents. 

We examine two separate aspects of this problem. 
Firstly, we consider the points of the phase diagram 
shown in Fig. [21 for a fixed population of the two gases. 
Secondly, we fix the interaction strengths and vary the 
ratio N A /N B . In both cases we examine the imaginary- 
time evolution of initial states with some finite, nonzero 
expectation value of the angular momentum, in the ab- 
sence of any external rotation of the trap. For a given 
set of parameters, the existence of a converged final state 
with a nonzero expectation value of / implies that the 
associated current is metastable. 



1. Variable couplings and fixed populations 

We have examined all the points that are shown in 
the phase diagram of Fig. [5] Those corresponding to a 
final state with a nonzero expectation value of the an- 
gular momentum are represented in Fig. 2 with "solid" 
symbols, while the ones that decay to the non-rotating 
ground state are represented with "open" symbols. 

Clearly, the vast majority of states correspond to 
metastable currents, except those close to the diagonal 
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9 = 9ABi as we ll as those with sufficiently small g (for 
all values of qab)- The obtained results show that the 
angular momentum of the metastable states is always an 
integer multiple of the particle population, which implies 
that the associated densities are necessarily circularly 
symmetric. This is consistent with the statement that 
circular symmetry is a necessary - though not sufficient 
- condition for the (meta)stability of the currents, as oth- 
erwise the circulation may escape from the gas (since in 
this case there is no barrie r sep arating the rotating state 
from the non-rotating one [181]). 



2. Variable relative population and fixed couplings 

Let us now study the effect of a variable relative pop- 
ulation between the two gases on the (meta) stability of 
the currents. We thus fix the interaction strengths, as 
well as the population of the one component {Nb), and 
study this question starting from the case with Na = 0, 
all the way up to the limit Na 3> Nb, with an initial 
state that has some angular momentum in component 
B. Since Nb is fixed, the above procedure corresponds 
physically to keeping the scattering lengths fixed. 

As mentioned above, the (meta)stability of the per- 
sistent currents depends on the competition between az- 
imuthal phase separation and circular symmetry of the 
density distribution of the two components. Thus, when 
Na = 0, the B component spreads within the whole sys- 
tem for any value of g and the currents can be metastable, 
provided that the coupling g is sufficiently large. If the 
population of the component A becomes nonzero but is 
small enough, it acts only as a weak perturbation, and 
azimuthal symmetry is still preserved. 

However, beyond a critical ratio Na/Nb, the two 
species separate azimuthally and the currents are no 
longer metastable. This critical value depends on the 
actual intra- and inter-component interaction strengths. 
A further increase of Na with respect to Nb drives the 
system to a phase of radial separation, and metastability 
is recovered. 

Finally, in the limit where Na becomes too large, this 
component cannot fit into only one of the two rings. This 
leads again to azimuthal phase separation of the gases, 
which is preserved even in the limit Na 3> Nb- As a 
result, the stability of the currents is lost again. 

We show in Fig.[S]the densities for the case g = 35, and 
9ab = 55 and for various values of the ratio Na/Nb, 
in order to illustrate the mentioned sequence of phase 
transitions. In particular, the results correspond to 
Na/Nb = 0.01, N A /N B = 0.25, N A /N B = 0.75, and 
Na/Nb = 10, with the first and the third cases corre- 
sponding to the metastable states. According to our sim- 
ulations, changing the values of the interaction strengths 
modifies the sizes of the "windows" in Na that separate 
the different phases, but yields qualitatively similar re- 
sults. 




FIG. 5: (Color online) Density distributions ua{p,0) = 
\<j>A(p, 0)\ 2 and n B (p,9) = \<f> B (p,0)\ 2 with g = 35, and 
gAB = 55 of the component A (dark blue) and B (light 
yellow), for (a) N A /N B = 0.01, (b) N A /N B = 0.25, (c) 
Na/Nb = 0.75, and (d) Na/Nb = 10. In the cases (a) and 

(c) the states support metastable currents, while in (b) and 

(d) they do not. The peaks in the densities correspond to the 
minima of the inner and the outer parabolae of the confining 
potential. 



V. SUMMARY AND CONCLUSIONS 

As shown in the present study, a coupled system of 
two distinguishable Bose gases that interact with an 
effectively-repulsive contact potential and are loaded in 
a concentric double annular trap reveals a series of phase 
transitions and the existence of metastable currents. 

For weak interactions, when the chemical potential is 
much smaller than the barrier that separates the two po- 
tential minima, the gases are confined in the inner ring, 
with the width of their transverse profile being smaller 
than the radius of the ring. Thus, their motion is (at 
least) close to being quasi-one-dimensional. On the other 
hand, as the couplings increase, the barrier plays a de- 
creasingly important role, and their transverse width be- 
comes comparable to the radius of the ring(s). Such a 
trapping potential interpolates up to some extent be- 
tween one-dimensional and two-dimensional motion, de- 
pending on the strength of the coupling between the 
atoms. The interplay between this effect and the strength 
of the inter- and intra-atomic couplings gives rise to inter- 
esting phase transitions in the ground state of the system, 
including three different geometries: phase coexistence, 
radial phase separation, and azimuthal phase separation. 

An interesting feature of the system considered is its 
response to some finite rotational frequency of the trap. 
The basic picture resembles very much the one where the 
couplings are varied, inducing axial and/or radial phase 
separation of the two components. 

The robustness in the (meta) stability of the currents 
that we found for the vast majority of the points in the 
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phase diagram of Fig. [2] is another interesting aspect of 
this study. Metastability is not found for the cases when 
the couplings are all nearly equal or exactly equal to each 
other, as well as for the cases with small g, independently 
oig A B- 

According to Ref. 20], a necessary (but not sufficient) 
condition for metastability is that the trapping potential 
does not increase monotonically from the center of the 
trap. The present results suggest that multiple variations 
in the monotonicity of the trapping potential enhance the 
stability of the currents. 

Last but not least, the discontinuous phase transitions 
we have found both in the ground-state phase diagram 
when the two gases separate radially as the couplings are 



varied, as well as in the response of the system when 
the rotation of the trap is varied, imply that hysteresis 
should show up as the coupling/rotational frequency in- 
creases /decreases. 
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